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ABSTRACT 

The origin of R Coronae Borealis stars has been elusive for over 200 years. Currently, two theories for their formation have been 
presented. These are the Final Flash scenario, in which a dying asymptotic giant branch (AGB) star throws off its atmosphere to reveal 
the hydrogen poor, heavily processed material underneath, and the double degenerate scenario, in which two white dwarfs merge 
to produce a new star with renewed vigour. Some theories predict that the temperatures reached during the latter scenario would 
destroy any lithium originally present in the white dwarfs. The observed lithium content of some R Coronae Borealis stars, therefore, 
is often interpreted as an indication that the Final Flash scenario best describes their formation. In this paper, it is shown that lithium 
production can, indeed, occur in the merging of a helium white dwarf with a carbon-oxygen white dwarf if their chemical composition, 
particularly that of 3 He, is fully considered. The production mechanism is described in detail, and the sensitivity of lithium production 
to the merger environment is investigated. Nucleosynthesis post-processing calculations of smoothed-particle hydrodynamics (SPH) 
tracer particles are performed to show that any lithium produced in these environments will be concentrated towards the cloud of 
material surrounding the R CrB star. Measurements of the lithium content of these clouds would, therefore, provide a valuable insight 
into the formation mechanism of R CrB stars. 
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1. Introduction 

Some chemically peculiar stars are deficient in hydrogen, while 
being enriched in carbon and oxygen. They are known col- 
lectively as hydrogen-deficient stars, and can be broken into 
three sub-groups: Hydrogen Deficient Carbon (HdC), Extreme 
Helium (EHe), and R Coronae Borealis (R CrB) stars. Their 
origin is difficult to explain through standard stellar evolution 
theory, so alternative scenarios for their formation are needed. 
One such theory is the Final Flash (FF) scenario, in which 
a late helium-shell flash occurs in a post-Asymptotic Giant 
Branch (AGB) star, hence moving it back towards the AGB part 
of the Hertzsprung-Russell diagram. The numerical models of 



Schonberner ( 1979b first showed that fully developed helium 



shell flashes at temperatures greater than 100 MK are possi- 
ble, lending weight to this theory. Indeed, FF objects, such as 
Sakurai's object, are known to exist (Duer beck & Benetti|1 996). 
The other leading theory is known as the double degenerate (DD) 
scenario, which involves the merging of two white dwarf stars - 
one with a carbon-oxygen (CO) core, the other, helium (He). Full 
understanding of these two production mechanisms is yet to be 
reached, and considerable effort is still needed to explain all of 
the observational signatures of these stars. 

A valid question that arises when considering the DD sce- 
nario is as follows: How feasible is it that two white dwarfs 
with the required compositions and masses merge to form an 
object resembling an R CrB star? To answer this question, the 
rate of these events within an observable distance can first be 



estimated. Roughly two of every three stars are born in binary 
systems ( Nelemans et al.||2001| l. Of these stars, however, only 
a fraction are expected to evolve into double white dwarf sys- 
tems that will merge within one Hubble time. Various estimates 
for this fraction are available, but as an example, the models of 
|Nelemans et al. ( 2001) > suggest that the birth rate of double white 
dwarfs in our galaxy is around 0.05 yr _1 , of which around 50% 
will merge within one Hubble time. About 15% of those merging 
systems will comprise of a carbon-oxygen and a helium white 
dwarf (where the CO white dwarf was formed first). While direct 
comparisons of this rate with the population of known hydrogen 
deficient stars is not trivial, these rates are consistent with the 
known population of close binary systems ( jlben et al.|1997| l. 

The DD scenario is also successful, at least qualitatively, in 
explaining the bulk of the surface abundances in hydrogen defi- 
cient stars (e.g., |Jeffery et al. |201 1 1. Recently, models with de- 
tailed nuclear networks have shown that nucleosynthesis can oc- 
cur in the DD scenario, thus explaining the over-ab undance of 



19 



F and O in R CrB stars (Longland et al. 



2011 
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However, 

some details such as the enrichment of s-process elements are 
more difficult to explain. Lithium abundances also pose a prob- 
lem, with four stars analysed by Aspl und et al.| ( |2000[ ) yielding 
over-abundances in the range of 1.4 to 2.3 dex. It has been hy- 
pothesised that lithium can be enriched to these levels success- 
fully in the FF scenario, while being hard to reconcile with high 
oxygen abundances in the DD scenario (Aspl und et al.| [2000; 
Clay ton et al.|201l) . The authors of those studies therefore in- 
terpreted the success of lithium production in the FF scenario as 
evidence that hydrogen deficient stars are most likely created in 
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the dying stages of AGB stars. However, the formation of lithium 
in the DD scenario has not been fully explored to date. 

The aim of this paper is to show that lithium can, indeed, be 
produced in the merging of a helium white dwarf with a carbon- 
oxygen white dwarf. We will begin in Sect. [2]by discussing the 
evolution of main sequence stars into a white dwarf binary sys- 
tem before explaining the mechanism necessary for lithium en- 
richment in Sect. [3] A number of numerical tests of this produc- 
tion are presented in Sect. |4] and its sensitivity to the thermody- 
namic properties of the merging event is discussed. These nu- 
merical tests will be compared with a detailed nucleosynthesis 
study by using Smoothed Particle Hydrodynamics tracer parti- 
cles in Sect. |5] and we will conclude in Sect.|6] 



2. Prior evolution of the binary system 

In order to understand the mechanism responsible for producing 
lithium in the DD scenario, the evolution of the binary system 
must first be addressed. The most likely system giving rise to 
the necessary white dwarf pair (one 0.6 M carbon-oxygen and 
one 0.4 M helium white dwarf) is a system containing stars with 
masses of 5 M Q and 1 M , respectively. Following two common 
envelope episodes, mass is ejected from the system resulting in 
two white dwarf stars in a close orbit (see |Iben & Tutu kov 1984; 
|Nelemans et al.|2001| for more information). 

Before considering the nucleosynthesis that occurs during 
the merger, the compositions of the white dwarfs must be un- 
derstood. It has been pointed out (Saio & Jeffery 2002; Jeffery 
etal.|20lT) that the surface composition of merging white dwarfs 



is different from their core compositions. This is because during 
the mass loss phases of the stars, not all unburned material is 
lost in stellar winds, thus thin, partially burned layers remain on 
the white dwarf surfaces. Therefore, some buffer regions were 
defined in those studies to take this material into account. This 
material is partially processed because it was originally close to 
the hydrogen burning shell in the progenitor star, so its chemical 
make-up should be altered from its original composition. 

To estimate the chemical abundances in the partially pro- 
cessed layer, let us consider how hydrogen burns in these en- 
vironments. For the helium white dwarf (originally a 1 M Q star), 
hydrogen should be consumed through the CNO cycle (dominat- 
ing the energy production) and the pp-chain, which converts four 
protons into helium while releasing energy in the form of pho- 
tons and neutrinos. The pp-chain follows three paths depending 
on the temperature of the environment. The dominant proton- 
buming paths with respect to ascending temperature are: 

PPI p(p,jS + )d(p,y) 3 He( 3 He,2p) 4 He, 

PPII p(p,jS + )d(p,y) 3 He( 4 He,y) 7 Be(e-, v) 7 Li(p, 4 He) 4 He, 

PPIII p(p,j6 + )d(p,y) 3 He( 4 He, r ) 7 Be(p,7) 8 B(y3 + ) 8 Be^2 4 He. 

A series of partial differential equations can be written to de- 
scribe these chains (see, for example, Iliadis 2007). Of particu- 
lar interest to us are those governing the evolution of 3 He as will 
become apparent in Sect. [3] By assuming that deuterium quickly 
reaches its equilibrium abundance in the star, the equation gov- 
erning 



3 He evolution can be written as 



d( 3 He) (H) 2 



dt 



<OT>p 



( 3 He)Vv> 33 



( 3 He)( 4 He)<<rv) 34 , (1) 



where (H), ( 3 He), and ( 4 He) are the abundances of hydrogen, 
3 He, and 4 He, respectively; (crv) pp is the rate of the p(p,/? + )d 
reaction; (<xv)33 is the rate of the 3 He( 3 He,2p) 4 He reaction; and 
(<xv)34 is the reaction rate of 3 He( 4 He,y) 7 Be. The equilibrium 



abundance of 3 He, ( 3 He) e , is determined by finding the roots of 
this quadratic equation when <f( 3 He)/df = 0, the only physical 
solution being 



( 3 He) f = — — [-( 4 He)< C rv>34+ 

2<crv>33 L 

<7v)pp{<rvh-$ + ( 4 He) 2 <<xv) 



(2) 



Hydrogen shell burning in a red giant branch star with approx- 
imately 1 M G occurs at 10-20 MK, thus energy is produced pri- 
marily through the CNO cycle. However, the pp-chain is still 
active in these conditions with a minor contribution to energy 
production and is essential for the production of 3 He, and hence 
7 Li. Throughout the following discussion, a point in the hydro- 
gen burning shell defined as the position at which H = 0.5 and 
( 4 He) = 0.5 is considerecQ This point, with an estimated temper- 
ature of T — 15 MK, is destined to become part of the hydrogen 
buffer in both white dwarfs once their envelopes have been lost 
to stellar winds. By using these values in Eq. O), 3 He mass frac- 
tions of ( 3 He) e = 1 x 10~ 5 are found to be reached in 5x 10 5 years, 
well below the red giant branch lifetime for these stars. Note 
that the bulk of 3 He is produced in the envelope of the star at 
much larger radii (and lower temperatures) than considered here. 
At later stages of the star's evolution, the envelope material is 
convectively mixed into this region, consequently increasing the 
3 He content relative to what is produced in the hydrogen-burning 
shell. In this section, we do not take these time-dependent effects 
into account, but rather take the conservative approximation that 
no convective mixing takes place to increase 3 He in the partially 
H-burned shell. These approximations are only made in Sects. [2|- 
|4]in order to illustrate the mechanism for producing 7 Li in white 
dwarf mergers. In Sect. |5j these approximations are relaxed to 
take the full evolutionary history of the white dwarfs into ac- 
count. In the centre of the hydrogen burning shell in a 1 M red 
giant star, one would therefore expect 3 He to be overabundant. 
Later, once the star's envelope has been eroded during the com- 
mon envelope stage, this 3 He persists in the surface layers of 
both white dwarfs. The fate of this material becomes essential in 
understanding lithium production during white dwarf mergers. 



3. Merging event 

Following the common envelope stages, gravitational wave 
emission is responsible for further reducing angular momen- 
tum in the system until finally, once tidal forces become strong 
enough, the helium white dwarf is disrupted rapidly (in about 
100 seconds) and is accreted on to the carbon-oxygen white 
dwarf. This merging episode results in an object that consists 
of three parts: (i) a dense, central object consisting of the core of 
the carbon-oxygen white dwarf; (ii) a hot corona in which nu- 
cleosynthesis has occurred; and (iii) a thick accretion disk (e.g., 
|Loren-Aguilar et aL 2010) . While the hot corona should con- 
tain some material originating in the outer regions of the carbon- 
oxygen white dwarf, thus explaining the high abundances of car- 
bon and oxygen in hydrogen deficient stars, it consists primarily 
of material from the helium white dwarf. 

As we have shown in Sect. [2j the hydrogen rich buffer of the 
white dwarfs is expected to be enriched in 3 He. This is of par- 
ticular interest because helium burning conditions are reached 



1 Note that in reality, the outer surface of the helium white dwarf does 
not consist of a discontinuity of abundances, but rather a continuous 
profile of abundances with respect to radius 
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during the merger, so the unstable nucleus 7 Be can be synthe- 
sised, which subsequently captures an electron to produce 7 Li. 
This synthesis depends on rapid heating and cooling of mate- 
rial on timescales that are rare in astrophysical environments 
(see Sect. [4]). The production reaction chain itself is simple: 
3 He(a,y) 7 Be(e~, v) 7 Li. It is complicated somewhat because 7 Be 
can be destroyed by proton capture to produce 8 B, which sub- 
sequently decays back into two a-particles or undergo further 
proton capture to produce 9 C. Furthermore, at high tempera- 
tures, the reverse reacti on 8 B(y,p) 7 Be is particularly importan t 



Table 1. Sensitivity study variables 



for 7 Be production (e.g.,|Boffin et al. 



1993 



Hemanz et al 



1996) 



equa- 



Nevertheless, at low temperatures, the partial differentia 
tions discussed in Sect. [2 are consulted once again to provide us 
with the 7 Be equilibrium abundance: 



( 7 Be) e 



( 3 He) e ( 4 He)<crv) 3 4 
A~l + H{crv)\i 



(3) 



where Aj is the electron capture rate of 7 Be into 7 Li and (crv)n 
is the rate of the 7 Be(p,y) 8 B reaction. 

The behaviour of this equilibrium abundance and its 
timescale are important in 7 Li production. As the stellar temper- 
ature rises, so does the ( 7 Be) c abundance, and the time required 
to reach those equilibrium values decreases dramatically. At a 
temperature of T = 100 MK, therefore, equilibrium abundances 
of 7 Be e = 1 x 10~ 7 are reached in about one second for material 
originating from the centre of the hydrogen buffer. Therefore, if 
high temperatures are reached followed by rapid cooling before 
exhaustion of the 3 He reservoir (the time-scales involved do not 
allow for replenishment of 3 He), 7 Be will freeze-out at highly 
enriched values. This 7 Be is then free to undergo el ectron cap- 
ture, p roducing 7 Li with a half-life of about 53 days (Audi et al. 
2003} , resulting in a stellar environment enriched in Li. 

We have shown that, in principle, 7 Be (and consequently 
7 Li) can be produced in high concentrations if the environment 
reaches the required temperatures for the correct time period. 
The obvious question that arises, therefore, is: how sensitive 
is this 7 Be production to the thermodynamic properties of the 
merging material? 



4. Numerical tests 

In the case of a white dwarf merger, the buffer layer of the he- 
lium white dwarf is expected to be amongst of the first material 
to impact the surface of the carbon-oxygen white dwarf, thus 
reaching high temperatures - typically above 300 MK ( |Loren- 
Aguil ar et al.p OlO). The energy generated by nuclear burning 
at these temperatures ejects the processed material back out to 
larger orbits, cooling it and allowing for any 7 Be to decay into 
7 Li. The assumed temperature density profile for buffer material 
falling on to the surface of a carbon-oxygen white dwarf there- 
fore consists of three stages: (i) a rapid exponential rise; (ii) a 
brief stage in which both the temperature and density remain 
constant; and (iii) an exponential, adiabatic cooling stage. 

A small hydrogen burning network is used consisting of nu- 
clei from A=l to A=12 to study the 7 Li production sensitivity on 
the profile parameters. The network contains all reactions perti- 
nent to nucleosynthesis in the pp-chain, including reverse reac- 



tions. Rates are adopted from the REACLIB database (Cyburt 
et al.pOlO i. In order to approximate the production of 'Li in 
the merger, the same location in the hydrogen buffer as used in 
Sect.[3]is considered. This region has 'H and 4 He abundances of 
(H) = ( 4 He) = 0.5 and a 3 He abundance of ( 3 He) = 1 x 10~ 5 , 



Variable 


Minimum value 


Maximum value 


tr(s) 


0.01 


20 


tp (s) 


0.01 


50 


to (S) 


0.1 


50 


Tmax (MK) 


80 


800 



Notes. Summary of the variable ranges used for the 7 Li production sen- 
sitivity study. refers to the rise time to maximum temperature, T ma x- 
t p is the time spent at this peak temperature, and t c is the time taken to 
cool back to T = 1 MK. 



as calculated in Sect. [2] All other nuclei have negligible initial 
abundances. 

The 7 Li production sensitivity is investigated with respect to 
four profile parameters: (i) rise time, t r , taken to heat the mate- 
rial from T = 1 MK to the maximum temperature, Tmax; (ii) 
cooling time, t c (i.e., the time taken to cool to T — 1 MK); 
(iii) amount of time spent at peak temperature, t p ; and (iv) max- 
imum temperature Tmax- 7 Li production is also sensitive to the 
initial abundances, depending mostly on the initial concentra- 
tion of 3 He. The parameter ranges adopted are summarised in 
Table [T] The density is chosen to vary from po = 1 x 10 3 g/cm 3 
to pmax = 1 X 10 5 g/cm 3 , with a cooling tail that follows an 
adiabatic cooling curve. This behaviour is consistent with the 
SPH models of Loren-Aguilar et al. ( 2010) . An equally spaced 
grid of 20 values are used for each variable, thus nucleosynthe- 
sis arising from 160000 profiles is analysed. A maximum 7 Li 
abundance of ( 7 Li) = 1.4 X 10~ 5 (corresponding to a factor of 
1000 more lithium than is found in our solar system) occurs at 
temperatures of 350 MK with the fastest rise, cooling, and peak 
times considered. These extreme values are, however, unrealis- 
tic. In reality, rapid heating of the material is expected, but the 
cooling time is subject to variability depending on the dynamical 
situation of the merger. We therefore consider the dependence of 
7 Li production to peak temperature and cooling time assuming 
a nominal rise and peak times of 1 s. This sensitivity is shown 
in Fig. [T] A number of interesting features are apparent. First, 
lithium is produced most efficiently for sharp profiles (i.e., short 
rise and cooling times). This is because of the freeze-out mech- 
anism discussed in Sect. [5] coupled to the fact that 3 He replen- 
ishment cannot occur quickly enough to maintain the 3 He: 7 Be 
equilibrium at peak temperatures. Secondly, the production of 
lithium as a function of temperature seems to follow a rather 
complicated pattern that requires some detailed discussion. 

By considering Fig. [T] three temperature regions are defined 
according to how much Ti is produced: (i) the low temperature 
region in which 7 Li production is independent of fall time; (ii) 
the medium temperature region of 200 < Tmax < 400 MK, in 
which the 7 Li abundance depends strongly on the profile shape 
and peaks around Tmax = 300 MK; and (iii) the high tempera- 
ture region in which very little 7 Li is produced. 

In the low temperature regime, around Tmax = 100 MK, 
the rate of the 3 He(a-,y) 7 Be reaction is not high enough to fully 
destroy 3 He during the merger. Therefore, an equilibrium is 
reached, which is rather insensitive to the dynamics affecting 
the shape of the profile. The equilibrium 7 Be abundance then 
follows the behaviour discussed in Sect. [3] However, as the tem- 
perature rises, the reaction rates increase, thus causing 3 He to 
be partially, and eventually fully destroyed. In this event, 7 Be 
can only be replenished by the photo-disintegration reaction, 
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Fig. 1. (Colour online) Sensitivity of 7 Li production to the cool- 
ing time, t c , and maximum temperature, Trnax, of the profile. 
The contour labels refer to log lo ( 7 Li)/( 7 Li) (i.e., zero repre- 
sents a lithium abundance equal to solar values). These values 
have been calculated assuming a rise time of t r — 1 s and a peak- 
time of t p — 1 s. 



8 B(y,p) 7 Be, whose rate is low around 200 MK, thus explaining 
the rapid drop in 7 Li production from 150 MK to 200 MK. 

In the medium temperature regime, 7 Be can no longer be 
produced from 3 He, and the equilibrium ratio of 7 Be/ 8 B is ap- 
proximately proportional to the ratio of their reaction rates: 
^i/io'vhs (where /t S7 is the photo-disintegration rate of 8 B 
into 7 Be and (crv)^ is the proton capture rate of 7 Be(p,y) 8 B). 
Assuming that all of the 3 He is converted to 7 Be and 8 B early 
in the merging event, and noticing that the reaction rate ratio 
increases with temperature, it follows that the surviving 7 Be 
should increase with temperature. However, one further effect 
needs addressing: the destruction of 8 B by the 8 B(p,y) 9 C reac- 
tion. This reaction, active above T = 100 MK, erodes the 8 B 
abundance and thus the abundance of 7 Be. For this reason, longer 
profiles (i.e., profiles that take longer to cool) suffer from lower 
7 Li production. 

At high temperatures, above T = 400 MK, another process 
comes into play: the destruction of 7 Be by the 7 Be(o',y) 11 C reac- 
tion. This reaction serves to divert the reaction flow away from 
8 B, and thus, the equilibrium 7 Be/ 8 B cannot be established. The 
rapid increase of the 7 Be(a,y) n C rate as a function of tempera- 
ture causes the sharp cut-off in 7 Li production shown in Fig. [I] 
Therefore, in helium rich environments, 7 Be is not expected to 
survive destruction by this mechanism at temperatures above 
T = 400 MK. 



5. Results from hydrodynamic models 

5. 1 . White dwarf models 

As shown in the preceding sections, the initial composition of the 
white dwarfs is critical to understanding lithium production dur- 



ing the merger event. Therefore, the detailed white dwarf mod- 
els of |Althaus et al.| ( |2010) and |Renedo et al.| ( |20T0) are used to 
evolve stars from the main sequence to the white dwarf stage. 
These evolutionary models contain a limited nuclear network of 
16 elements: 'H, 2 H, 3 He, 4 He, 7 Be, 12 C, 13 C, 14 N, I5 N, 16 0, 17 0, 
ls O, 19 F, 20 Ne, and 22 Ne. This network is sufficient for following 
nuclear energy generation during the evolution of the star. 

The carbon-oxygen white dwarf model is obtained by evolv- 
ing 2.5 M Q main sequence stars with metallicities of Z = 1.5 x 
10~ 2 and Z = 1 x 10~ 5 through the thermally pulsing AGB 
stage and subsequent mass loss period (using the prescription 
of Vassiliadis & Wood 1993 i through the planetary nebula phase 



to the white dwarf phase. The subsequent cooling of the white 

for more 



dwarf is carefully modelled (see Althaus et al 



2010 



information) to take into account the effects of, among others, 
diffusion of material. 

The helium white dwarf model is computed similarly to the 
carbon-oxygen white dwarf. However, since the orbital radius of 
the two stars decreases during the first common envelope stage 
(i.e., when the more massive star loses mass in its AGB phase), 
the common envelope could occur earlier in the stars evolution- 
ary path. A 2.5 M main sequence star is therefore used. During 
the stars evolution on the red giant branch, an artificially high 
mass-loss rate is applied to account for the common envelope 
stage. As a consequence of this high mass-loss rate, the star loses 
its envelope before core helium burning commences, and evolves 
on to the helium white dwarf cooling track. 

Following the evolution described above, white dwarf mod- 
els are obtained that describe the distribution of isotopes in the 
two stars as a function of radius. These models are then mapped 
onto our hydrodynamic models in order to follow the evolution 
of these abundances as the two white dwarfs undergo a merging 
event. 



5.2. Merging system models 

The hydrodynamic models used to provide us with input to our 
nucleosynthesis calculations are the same as those presented 
in |Guerrero et al.| p004|, ILoren-Aguilar et al.| ( |2009l [20To) i, 
and |Longland et aL ( |2011 i. A Lagrangian Smoothed Particle 
Hydrodynamics (SPH) code is used to obtain temperature- 
density profile input for the post-processing nuclear network. 
The model used for the calculations presented here is computed 
for the merging of a 0.4 M Q helium white dwarf with a carbon- 
oxygen white dwarf of mass 0.8 M Q . While this model does 
not correspond exactly with the white dwarf masses required to 
form an R CrB star, their dynamical merging nature should be 
similar. Additionally, since the envelope of the final system is 
mostly comprised of material originating in the helium white 
dwarf, this merging system should reproduce sufficiently the 
0.4 + 0.6 M Q case. The SPH particle masses were 2.6 x 10 -6 M Q 
and 5.3 x 10~ 6 M o for the two stars, respectively to provide a 
total of 300 000 particles. 

A sample of 10000 tracer particles is chosen to follow the 
chemical evolution of material during the merger. Following the 
merging event, only material that is present in the hot corona 
and surrounding accretion disk will be visible to any observer. 
Consequently, only particles that finish their evolution with or- 
bital radii outside the dense central object are chosen as tracers. 
These particles correspond to those that would form the atmo- 
sphere and chromosphere (Clayton 1996) of a hydrogen deficient 
star. A similar nuclear network as that presented in Sect.|4](with 
additional isotopes introduced to account for those in our white 
dwarf models) is used, which is sufficient to model the synthesis 



4 



Richard Longland et al.: Lithium production in the merging of white dwarf stars 










1 


- 


Hot 


1 Accretion 






Corona 


I Disk 






y\ It \ 
\ r \ 

y ' * N 1 \ 


J 





4e-03 



4e-02 



4e-01 



R R. 



Fig. 2. (Colour online) Lithium abundance yield obtained from 
our SPH models as a function of radius in the post-merger ob- 
ject. The abundances are represented as logarithmic abundances 
relative to solar values i.e., [Li] = log 10 ((Li)/(Li) G ) and (Li) is 
the mass fraction of lithium. R/R represents the radius in units 
of solar radii. The solid black line represents the radial lithium 
abundance resulting from the merging of solar metallicity stars, 
while the dashed red line represents that of Z = 1 x 10~ 5 stars. 
The single point with error bars represents the range of observed 
lithium abundances in R CrB stars (see |Jeffery et al.||201 f and 
references therein). 



of light elements. Elemental abundances, for comparison with 
observational data, are then calculated by summing the isotopic 
abundances obtained from the models. 



5.3. Model results 

Average lithium abundances, normalised to solar values, as a 
function of radius in the merger product are shown in Fig. 2] 
Clearly, lithium enrichment is present in the outer regions of the 
object with respect to solar abundances. The radii at which this 
high enrichment are found are, indeed, rather large, and corre- 
spond most likely to the accretion disk surrounding the merger 
product. The reason for this enrichment is that the buffer ma- 
terial, which is destined to produce lithium because of its high 
concentration of 3 He, is amongst the first matter accreted on to 
the carbon-oxygen white dwarf. Consequently, it will heat up 
rapidly and is forced out to larger orbital radii. Subsequent ac- 
cretion of helium rich matter from the destroyed helium white 
dwarf prevents this material from falling back on to the carbon- 
oxygen white dwarf's surface and it remains in the accretion disk 
following the merging event. 

Lithium production is also expected to be independent of 
metallicity (it depends mainly on the 3 He enrichment of the 
white dwarfs). Indeed, this is the case in the low metallicity mod- 
els considered (represented by a dashed red line in Fig. [2j where 
the lithium abundance found in the outer regions of the final ob- 
ject agree very well with those found with our high metallicity 
models. 

The distribution of lithium in our models could explain why 
lithium is only observed in a limited number of R CrB stars. 
Although no accretion disks have been detected, our SPH mod- 
els suggest that a rapidly rotating Keplerian disk should be 



formed subsequent to the merging event (see Loren-Aguilar et al. 
2010, and references therein). This disk could be partially re- 
sponsible for the observed circumstellar material surrounding R 
CrB stars. If this is the case, then if an observer were to view 
the object side-on, the accretion disk could obscure the atmo- 
sphere of the star and dominate any spectra obtained. If the star 
were face-on to the observer, however, radiation from the surface 
of the star itself will dominate the observers measurements, and 
hence, lithium abundance measurements would be rather low. 
Our discovery for the distribution of lithium in the merger prod- 
uct could also be used observationally to distinguish between the 
FF and DD merger scenarios. If, indeed, the DD scenario was the 
formation channel for the star, our models indicate that lithium 
will be present primarily in the cloud of material surrounding 
the star. The leading theory explaining the characteristic declines 
of these stars (Clayton 1996) suggests that the dramatic drop in 
luminosity of the star is caused by a "puff" of carbon-rich ma- 
terial that is being accelerated towards the observer. If the star 
can be observed early during one of these declines, the surface 
of the star will be partially eclipsed, a process known as "veil- 
ing" (Lambert et al. 1990), hence revealing the diffuse surround- 
ing material for measurement. If lithium is present in high quan- 
tities in this material, it could indicate that the DD scenario is, 
indeed, responsible for the formation of R CrB stars. 

6. Conclusions 

In this paper, we have investigated the production of lithium 
in hot mergers of helium and carbon-oxygen white dwarf 
stars. Both by performing analytical calculations and by post- 
processing nucleosynthesis of SPH tracer particles, significant 
lithium production is possible. 

In our analytical models, we have shown that this production 
is strongly dependent on the dynamical properties of the mate- 
rial in the merger and have placed constraints on the conditions 
required for lithium production. To produce large amounts of 
lithium the material, which was enriched in 3 He in prior evolu- 
tion, must reach high temperatures rapidly, and then cool rather 
rapidly also. Lithium production also depends sensitively on the 
maximum temperature reached by the material through a rather 
complicated interplay of reactions and their reverse counterparts. 
Maximum lithium yields are achieved at peak temperatures of 
around Tmax = 300 MK provided the material undergoes rapid 
heating and cooling episodes. 

Post-processing of tracer particles from SPH models of 
merging white dwarfs confirms that lithium can, indeed, be pro- 
duced in the DD scenario, thus satisfying the constraints investi- 
gated in our analytical models. By performing these calculations 
for the merging of a 0.4 M Q helium white dwarf with a 0.8 M G 
carbon-oxygen white dwarf, we have shown that the agreement 
between our models and the observed lithium surface abundance 
in hydrogen deficient stars is reassuring. Furthermore, our mod- 
els predict that lithium should be predominantly produced in the 
outer regions of the merger product. If measurements are made 
to identify the location of lithium in R CrB stars, they would 
provide a vital clue in determining the formation mechanism of 
these peculiar stars. Further studies should also be performed to 
investigate the influence of other factors such as the white dwarf 
masses, compositions, and dynamical properties on the nucle- 
osynthesis of lithium in white dwarf mergers. 

It was previously believed that lithium observations in some 
R CrB stars indicated that they are produced by a powerful final 
helium flash in a dying AGB star. As we have shown in this pa- 
per, lithium can be produced in the merging of two white dwarfs 
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if the conditions meet the criteria laid out here. Lithium cannot, 
therefore, be used as a strong indicator that the DD formation 
scenario is invalid for R CrB stars. 
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